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Abstract 

'  / 

( V 

A  brief  summary  of  work  completed  by  the  authors  during  the  period  9/1/79- 
8/31/80  on  the  investigation  of  dispersive  wave  mechanisms  for  enhancing  the 
conversion  efficiency  of  parametric  acoustic  arrays  is  presented  via  (1)  a  detailed 
review  of  the  physical  concepts  underlying  the  investigation,  (ii)  an  in-depth 
review  of  the  theoretical  models  upon  which  the  Investigation  la  based,  and 
(ill)  a  summary  of  results  obtained  to  date,  followed  by  a  statement  of  con¬ 
tinuing  objectives. 
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List  of  Symbols 

Rectilinear  Cartesian  Coordinates  and  time 
Velocity  Potential 
Excess  Pressure 

Circular  Frequency  (loq  being  an  arbitrary  reference  frequency) 

Sum  and  Difference  Frequency 

Velocity  Potential  Spectrum 

Excess  Pressure  Spectrum 

Density 

Low  and  High  Frequency  Speed  of  Sound  in  a  Monorelaxing 
Mediumi  respectively 

Honllnear  Coefficient  of  the  Equatlon-of-State 

Second-Order  Nonlinear  Coefficient  of  a  Fluid  (y*1+^  in  liquids) 

A 

Eigenmodes 

Wavenumber  at  Frequency  w 

s 

Transverse  Wavenumber 
Axial  Wavenumber 

Velocity  Potential  Spectral  Amplitudes 
Velocity  Potential  Weighting  Coefficients 

Fundamental  Mode  Primary  Wa.’e  Pressure  Amplitudes 

Net  Sum  or  Difference-Frequency  Phase  in  a  Dispersive 
Medium 

Attenuation  Coefficient  at  Frequency  o>s 

Sum  or  Difference-Frequency  Attenuation  Coefficient 

Net  Attenuation  Coefficient 


Mean  Primary-Wave  Attenuation  Coef  i  .lc  lent 
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Sum,  Difference,  or  Arbitrary  Half  Rayleigh  Distances, 
respectively 

Difference-Frequency  Directivity  Function 

Beam  Angle  relative  to  the.  direction  of  Wave  Propagation 

Projector  Area  normal  to  the  direction  of  Wave  Propagation 

Dispersivity  Parameter  in  a  Monorelaxing  Fluid 

Relaxation  time  in  a  Monorelaxing  Fluid 

Shear  and  Dilatational  Coefficients  of  Viscosity, 
respectively 

Coefficient  of  Thermal  Conductivity 

Specific  Heat  at  Constant  Pressure 

Thermo-Viscous  Absorption  Coefficient 

Low  and  High-Frequency  Wavenumbers  renpectively  in  a 
Monorelaxing  Fluid, 

Axial  Wavenumber  at  frequency  in  a  Monorelaxing  Fluid 

Modal  Pressure  Amplitudes  at  Frequency  w 

Pressure  Field  Weighting  Coefficients  in  a  Monorelaxing 
Fluid 

Gaussian  Beam  Spot  Size  at  Frequency  Ug 
Kronecker  Delta 

Particle  Velocity 
Mach  Number 

Rayleigh  Distance  at  Mear  Primary  frequency 
Scaled  Source  Level  Parameter 

Range 

Plane  wave  Scaled  Range  Parameter 


1.  Introduction 
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Parametric  Amplification,  as  originally  envisaged  by  Cullen^"  and  by  Tien 
2 

and  Suhl  ,  concerned  the  transfer  of  energy  to  a  weak  signal  of  frequency  d  via 
sinusoidal  perturbation  of  the  parameters  (e.g.  inductances  or  capacitances)  of 
an  electrical  transmission  line  at  frequency  2u).  In  effect,  if  a  weak  sig^a1  of 
frequency  u)  and  a  strong  pump  wave  of  frequency  2 to  are  simultaneously  present 
in  a  nonlinear  transmission  line,  interaction  occurs  between  them  giving  rise  to 
a  difference-frequency  component  (i.e.  2oi  -  u>  ■  u) which  augments  and  consequently 

3 

amplifies  the  signal  of  frequency  gj.  As  subsequently  shown  by  Roe  and  Boyd 
however,  the  presence  of  nonlinearly  generated  harmonics  greater  than  2m  significantly 
diminishes  and  ultimately  undermines  the  amplification  process.  For  example, 
nonlinear  generation  of  the  sum  frequency  (i.e.  2u>  +  u>  -  3m)  component  gives  rise 
via  degerative  interaction  with  the  second  harmonic  (i.e,  2m  -  3oj  ■  -m)  to  a 
depletion  of  the  amplification  gain  at  the  signal  frequency.  In  order  to  ensure 
the  efficacy  of  Parametric  Amplification  therefore,  filtering  circuits  muBt  be  applied 
to  the  transmission  line  to  make  the  system  appear  as  a  low  pass  filter  which  will  block 
frequencies  greater  than  the  second  harmonic.  In  the  case  of  nonlinear  electro¬ 
magnetic  transmission  lines  such  filtering  is  readily  achieved.  However, in  the  case 
of  nonlinear  optical  and  acoustical  parametric  wave  interactions  which  occur  in 
bulk  media  the  problem  of  minimizing  the  influence  of  degenerative  coupling  becomes 
extremely  difficult  to  realize.  For  this  reason  the  development  of  . ow  noise 
narrowband  Parametric  Acoustic  receivers  has  been  severely  inhibited.  Likewise, 
the  conversion  efficiency  of  Parametric  Acoustic  transmitters  (which  constitute 
a  generalization  of  the  basic  form  of  Parametric  Amplification  described  above) 

A 

as  envisaged  by  Westervelt,  where  a  difference-frequency  signal  m  Is  formed  via 
nonlinear  interaction  of  high  intensity,  high  frequency  pump  waves  of  frequencies 
and  0)2  In  a  bulk  medium,  is  severely  diminished  by  degenerative  coupling  effects. 

In  this  instance,  the  conversion  efficiency  is  reduced  both  by  pump  depletion 
via  energy  transfer  to  nonlinearly  generated  higher  harmonics  (i.e.  2nu>^,  Znu^) 


5. 


and  by  degenerative  coupling  between  these  harmonics  and  the  upper  sideband 
int.ermodulation  frequency  components  (i.e.  noi^  +  mfi^)  •  As  shown  by  Tjotta^,  the 
most  significant  of  the  latter  interactions  is  that  which  occurs  between  the  com¬ 
paratively  strong  sum  -  frequency  component  and  the  secondary  pump  wave  harmonic 

One  solution  to  the  problem  of  blocking  those  parts  of  a  nonlinearly  generated 
spectrum  that  inhibit  the  resonant  interaction  of  desired  spectral  components  in 
bulk  media  is  to  inhibit  the  amplification  of  unwanted  frequency  components  via 
dispersive  processes.  Since  each  spectral  component  travels  at  a  different  phase 
velocity  in  a  dispersive  medium,  the  amplitude  variation  of  spectral  components 
acquire  the  character  of  spatial  beats  due  to  the  accumulation  of  relative  phase 
shifts  during  the  course  of  propagation.  As  the  disperslvlty  increases  the  beat 
periods  decrease, thus  reducing  the  peak  amplitudes  of  the  spectral  components.  By 
ensuring  that  unwanted  components  of  a  nonlinearly  generated  spectrum  occur  in 
regions  of  strong  dispersivity  which  at  the  same  time  appears  virtually  dispe-sion- 
less  to  the  frequencies  of  interest,  only  the  latter  spectral  components  are 
strongly  coupled,  thus  ensuring  the  possibility  of  significant  parametric 
amplification.  On  account  of  the  strong  'inherent'  dispersivity  of  most  optical 
dielectrics  this  effect  has  been  successfully  exploited  In  nonlinear  optical 
parametric  amplifiers,  as  described  by  Bloembergen. ^  Since,  on  the  contrary, 
most  acoustical  media  are  weakly  dispersive,  very  few  cases  of  'dispersive  wave 

filtering'  have  been  realized  with  the  exception  of  Shiren's  '  beautiful  ex- 

aH~ 

perimental  induction  of  anomalous  dispersion  in  MgO  crystals  containing 

i|m|i 

of  f e  ions  via  applied  magnetic  fields.  Another,  exception  is  the  Parametric 

9  10 

Acoustic  Amplifier  realized  by  Ostrovskii  and  Papilova  ’  via  boundary-induced 

dispersion  in  fluid  filled  rectilinear  waveguides.  A  related,  but  unrealized  form 

of  boundary-induced  'dispersive  wave  filtering'  has  been  discussed  by 

11  12 

Zareobo,  Serdobol'skaya,  and  Chernobai  ’  .  In  this  instance,  plane  pump 
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waves  propagating  in  a  simple  resonator  of  length  L  are  reflected  from  frequency 

dependent  termination  impedances  such  that  upon  reflection  from  the  impedance 

at  L  the  resulting  sum-frequency  phase  shift  is  equal  to  tt  radians.  Since  the 

amplitude  of  the  sum.  frequency  component  is  therefore  zero  at  points  2L,  4L, 

etc.,  its  growth  and  resulting  degenerative  influence  will  be  significantly  reduced. 

Again,  the  use  of  inhomogeneties  (e.g.  bubbles)  to  effect  'dispersive  wave 

13 

filtering'  as  analysed  by  Zabolotskaya  and  Soluyan  appears  to  virtually  exhauBt 

the  extent  of  investigations  in  fluids  carried  out  to  date.  In  solids,  Parametric 

Amplification  based  on  the  interaction  of  pump  waves  propagating  at  oblique  angles 

with  respect  to  each  other  has  been  investigated  by  Zabolotskaya,  Soluyan,  and 

Khokhlov^"  ^ ,  Lord1^,  and  Ivanov  and  Pluzhnikov*"®,  based  on  earlier  theoretical 

work  .  As  described  by  Rudenko  and  Soluyan  ,  the  required  angle  for  'resonant 

interaction*  between  interacting  waves  is  determined  by  the  synchronism  conditions 

(i.e.  m  _  ■  nuj,  +  ,  k  ■  nk.  +  mk_).  For  waves  propagating  at  an  angle  these 

nm  l  —  z  _nm  _z 

conditions  are  not  satisfied  at  all  the  intermodulation  frequencies.  Since  longi¬ 
tudinal  and  transverse  waves  propagate  in  a  solid  with  different  velocities,  it 
is  passible  therefore  at  certain  intersection  angles  between  the  pump  waves  to 
satisfy  the  synchronism  conditions  at  the  difference  frequency.  In  such  instances, 
the  synchronism  conditions  are  violated  at  the  sum  frequency  and  hence  the  latter 
ia  effectively  suppressed.  This  in  turn  ensures  parametric  efficiency  enhancement 
via  reduced  degenerative  coupling  interactions. 

Summarizing  the  limitations  of  the  above  papers,  it  should  be  noted  that 
previous  investigations  of  parametric  gain  enhancement  via  'dispersive  wave  filter¬ 
ing'  have  been  peacemeal  and  generally  restricted  to  particular  cases  of  lossless 
plane  wave  propagation.  They  have  moreover,  been  based  on  the  ubiquituous  assumption 

that  the  effect  of  dispersivity  is  sufficient  to  reduce  all  nonlinear  coupling 

6  23 

to  that  of  three  frequency  (pump,  signal,  and  idler)  interaction  process.  * 

No  attempt  has  been  made  to  address  the  problem  of  how  much  dispersivity 


ia  required  in  the  presence  of  absorption  and  diffraction  losses  to  effect  varying 
degrees  of  parametric  amplification.  Of  course,  this  is  a  much  more  comprehensive 
and  difficult  issue  because  it  involves  the  analysis,  not  of  three,  but  of  all 
significant  nonlinear  wave  interactions  subject  to  absorption  and  diffraction 
losses,  which  must  inevitably  be  carried  out  via  numerical  methods.  It  is  our 
contention  however,  and  the  purpose  of  the  present  investigation  to  show  that 
no  meaningful  progress  can  be  made  in  utilizing  'dispersive  wave  filtering'  to 
enhance  the  conversion  efficiency  of  parametric  acoustic  arrays  until  this  matter 
is  resolved.  This  will  be  accomplished  by  computing  solutions  of  dispersive  non¬ 
linear  acoustic  wave  equations  in  terms  of  dimensionless  parameters.  The  secondary 
question  of  how  much  amplification  can  be  realized  in  particular  dispersive  media, 
in  particular  frequence  banda  etc.,  can  then  be  deduced  aa  a  consequence  of  these 
investigations. 


t 
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2.  Theory 


In  this  section  we  will  outline  the  methods  and  procedures  that  will  form 
the  basis  of  our  investigation,  including  some  new  analytical  results.  Although 
related,  we  wish  to  make  a  distinction  between  'boundary-induced*  dispersion  and 
'inherent'  dispersion,  the  former  being  self-evident,  and  the  latter  being  due 
to  relaxation  mechanisms. 

It  should  be  noted  that,  with  the  exception  of  Ostrovskii  and  Papilova' s10 

24-28 

work,  previous  investigations  of  finite-amplitude  wave  propagation  in 

bounded  media  have  not  been  concerned  with  the  problem  of  exploiting  dispersivity 

to  enhance  the  process  of  parametric  amplification.  Likewise,  in  the  case  of 

29-34 

inherently  dispersive  media  none  of  the  previous  work  has  been  concerned 

35 

with  this  question,  but  rather  with  the  process  of  soliton  formation  ,  which 
although  a  topic  of  great  interest  does  not  concern  us  here.  We  will  proceed 
therefore,  to  develop  an  analytical  procedure  that  links  the  two  dispersion  mechanisms 
under  consideration.  This  approach  is  similar  to  that  adopted  by  Bloejnbergen 
in  his  investigation  of  nonlinear  optical  wave  interactions  in  dispersive  media. 


2.1  Boundary- Induced  Dispersivitj 


In  order  to  investigate  the  process  of  enhancing  the  conversion  efficiency 
of  parametric  acoustic  interactions  via  boundary-induced  dispersivity  we  begin 
with  the  lossless  form  of  the  nonlinear  wave  equation  and  subsequently  introduce 
losses  in  the  frequency-domain  by  means  of  complex  wavenumbers.  As  given  by 


Blackstock  this  equation,  in  terms  of  the  velocity  potential  <$>,  assumes  the  form 


(V2  -  c~2  a2H  -  c“2  3t  (V(i>.Vcj>)  +  (Y  -  1)  Ot<t>)  cv2<h) 


) 

>  +  0 


c'2  5 

0  t 


•V<j>  +  +  0  (<|>3) 

2  c 

o  ' 


9. 


where  the  'substitution  corollary'  has  been  invoked  in  deducing  Eq.  (la)  froir 

Eq.  (1)  (i.e.  Ot<j>)  (V24>)  «*  c"2  (3t<p)  (3*4>)  -A  c'q2  Q^)2). 

For  the  case  of  wave  propagation  in  a  rectilinear  layer,  such  as  that 

depicted  in  Fig.  1,  if  e3  (x,y)  represents  the  m,n  linear  eigenraode  of  the 

tun 

structure  at  an  excitation  frequence  0)s  ,  and  is  the  corresponding  transverse 
wavenumber  for  this  mode,  then  by  definition 

(V2  +  (K™)2)  C,  -  0  ,  ?2  -  32  +  >1  (i 

where  (k”V  -  k2  -  (xf1)2  (2 
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In  this  notation  kgwU)g/co  and  xa  is  the  axial  wavenumber  corresponding  to  k^11  . 
We  now  express  $  in  terms  of  a  modal  expansion  of  the  form 
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where  the  unknown  coefficients  A3^  (z)  are  to  be  determined  by  substituting 
Eq.  (4)  into  Eq.  (la).  When  this  substitution  is  made  in  the  left-hand-side  of 


Eq.  (la)  we  obtain  via  Eqs.(2)  and  (3) 
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Once  the  eigenmodes  (x,y)  have  been  specified  for  particular  boundary  conditions 
S  ^  s*' 

the  B  coefficients  can  readily  be  evaluated, 

mnm  n.  ta  n 
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Since  p'  ■  -pq  3c<j)  +  0(<)i  ),  the  excess  pressure  spectrum  (x,y,z,)  is 


given  by  Eq.  (4)  as  p^  (x.y.z)  (x.y.z) 
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Assuming,  in  the  case  of  parametric  interaction  between  finite-amplitude  primary 
waves  of  frequencies  u)^  and  that  all  nonlinearly  generated  waves  other  than 
the  difference-frequency  (i.e.  -  u^)  signal  could  be  effectively 

suppressed  via  boundary-induced  dispersion,  Eq.  (10)  becomes 
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where  A  S  A^  ,  A  J  =  A  ‘  ,  and  a“‘  =  A®“  . 
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If  the  lowest  order  primary  wave  modes  alone  are  excited  and  A^,  A^  are 
both  constant,  which  in  this  instance  assumes  that  nonlinear  waveform  distortion 
is  small,  or  alternatively  that  the  primary  fields  are  only  subject  to  absorption 
losses  via  the  imaginary  parts  of  their  respective  wavenumber,  then  Eq.  (13c) 
becomes 
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Substituting  for  A^  in  Eq.  (4)  it  follows  that  the  difference-frequency  pressure 
field  is  given  by 
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In  order  to  teat  this  weak  finite-amplitude  solution  let  us  assume  that  the  boundary 
surfaces  of  the  layered  medium  are  removed,  but  that  the  medium  Itself  remains 
dispersive  due  to  some  inherent  physical  mechanism.  The  problem  now  corresponds 
to  the  case  of  plane  wave  parametric  interaction  in  an  inherently  dispersive  un¬ 
bound  medium.  Thus,  after  some  manipulation,  Eq.  (15)  becomes 
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where  is  the  nonlinear  coefficient  the  fluid. 
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Hence, Eq.  (16)  can  be  repressed  as 
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in  a  lossy  diapersionless  medium, (18c) 


where  k^,_  ■  -  k_  cos  8  •  2k_  sin4 (0/2),  and  oy*  2(*0  in  the  later  case. 

Eq.  (18b)  is  the  appropriate  result  for  a  weak  parametric  interaction  in  o 
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lossless  dispersive  medium,  being  similar  to  that  derived  by  Rudenko  and 
23 

Soluyan  for  the  second-harmonic  field  of  an  initially  monotonic  wave.  Like- 

iz 

wise,  if  the  field  described  by  Eq.  (18c)  is  multiplied  by  - —  D  (9)  ,  correspond- 

Z  (i)— 

ing  to  the  case  of  a  difference-frequency  signal  collimated  over  a  Rayleigh 

distance  zo  *  k_A0/2TT,  where  kQ  is  the  area  of  the  primary-wave  projector  and 

D  (6)  is  the  diffraction  pattern  of  the  aperture  at  the  difference-frequency, 

-a_z  -i2k_z  sin2(6/2) 

then  in  the  far-field  where  e  •  e  ■  0,  the  resulting  form 

37 

of  the  difference-frequency  signal  reduces,  as  required,  to  Maze  and  Tjotta's 

4 

modification  of  Westervelts'  solution. 

Although  the  solutions  of  Eqs.  (13a)  to  (13c)  considered  in  this  section 

are  only  approximate,  it  should  be  noted  that  exact  solutions  expressed  in  terms 

38 

of  elliptic  functions  exist  for  the  case  of  constant  coefficients  as  shown 
in  Appendix  A. 


2.2  Inherent  Diapersivity 

In  this  section  we  deal  with  the  case  of  parametric  array  formation  via 

nonllnearly  interacting  colinear  paraxial  waves  in  inherently  dispersive  media. 

If  such  an  interaction  occurs,  for  example,  in  a  monorelaxing,  thermo-viscous 

fluid,  the  governing  equation  for  the  velocity  potential  4>*  correct  to  second- 
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order  terms  is  given  by  the  following  modified  form  of  Eq.  (la): 
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-  1))  is  thermo-viscous  coefficient, 
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and  T  is  the  relaxation  time.  In  this 
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For  the  case  of  progressive  finite-amplitude  primary  waves  radiated  in  the 

redirection  by  a  single  projector  in  an  unbounded  medium,  previous  work  by 

40 

Woodsum  and  Westervelt  justifies  the  approximation 
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For  Che  case  of  a  square  aperture  radiator  located  in  the  plane  z  ■  0  centered  at 
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x  ■  0,  y  ■  0,  the  modes  e  are  given  by  Kogelnik  and  by  Cook  and  Arnoult  as 
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The  orthononnal  relation  satisfied  by  the  modes  is 


dxdy  e®  (x,y,z)  e®,  ,  (x,y,z)  -  <5  ,  6  ,  . 
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-*•6  (x-x')  <5(y-y'),  when  z"z'  (31s 

Assuming  again,  as  in  section  2.1,  for  the  case  of  a  weak  parametric  interaction 
between  finite-amplitude  primary  waveB  of  frequencies  U)^  and  w^ithat  all  non- 
linearly  generated  waves  other  than  the  difference-frequency  (i.e.  o)_  ■  -  u),,) 

signal  could  be  effectively  suppressed  via  Inherent  dispersion  in  the  medium, 

Eq.  (27)  becomes 
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where  ?  H  p  ,  P  11  =  P  ^  ,  and  P  '  a  P  1  ^  (3] 

mn  ran  mn  mn  mn  mn 

If,  as  in  section  2.1,  the  fundamental  primary  wave  modea  alone  are  excited  and 
U1 

both  P  ,  P  are  constant,  which  implies  that  nonlinear  waveform  distortion 
oo  oo  r 

is  small,  or  alternatively,  that  the  primary  waves  are  only  subject  to  viscous 
attenuation  and  spreading  losses,  then  Eq.  (32c)  becomes 
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from  Eq.  (28) 


Taking  the  Fourier  transform  of  Eq.  (22)  and  substituting  for  P^"  from  Eq.  (34) , 
the  difference-frequency  field  assumes  the  form 


U.y.z)  -  a' 


-L\J  L  P“-  (a)  eu-  <x*y*z> 

m,n»-“  mn  mn 


ib“?  p“o  p!o*  -4v'  f  “  .  *r, .  “VW 

2o  c4  v  ra,n— *  ran  ' 

2P0  ®  \lj  •'  0 


.W  UJ.  10* 

dx'dy'  e  “(x' ,y',z')  e  ^(x'.y'.z')  e  2  (x',y',z')  (35) 

mn  oo  oo 


From  5q .  (31)  this  becomes 


oi,  w  * 

r»  ^  n 


Prn  (x,y,z) 


iSu  P  x  P  ‘  -iv  z 

-  00  oo  AU)_ 


3 

2p  cJ 

'  O  ® 


f  ik 

dz' 


2tt(z-z'  ) 


u,  v 

dx'dy'  eq^(x' ,y' ,z’)  eo^(x' *y! .z')  e  ^  *  2 


18. 


But  Eq.  (29)  gives 
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To  the  beat  of  the  authors  knowledge  Eqs.  (41)  and  (42)  Is,  a  new  result  for 

the  propagation  of  a  paramatically  generated  wave  in  a  dispersive  medium.  If  this 

dlapersivity  Is  removed  and  k^  -  k^  *  k_  the  equation  reverts  to  the  form  invasti- 
44-47 

gated  by  Fenlon  in  hia  detailed  analysis  of  near  and  far-fleld  parametric 

acoustic  array  interactions.  As  it  stands  it  can  also  be  shown  that  in  the  far- 

fleld  Eqs.  (41)  and  (42)  reveals  a  shift  In  the  angular  spectrum  dependent  on  the 

30 

amount  of  dlapersivity  which  was  first  referred  to  by  Novikov  and  subsequently 
48 

by  Fenlon.  Without  going  into  further  details  now  however,  it  is  our  intention 
to  fully  evaluate  Eqs.  (41)  and  (42)  numerically  in  order  to  determine  the  nature 
of  the  spreading  losses  so  that  the  complete  set  of  spectral  equations  may  be 
solved  along  the  beam  axis.  In  this  manner  the  effect  of  dlapersivity  in  destroy¬ 
ing  unwanted  spectral  components  will  be  investigated. 

Finally,  on  other  task  that  must  be  carried  out  is  to  elaborate  on  the  form¬ 
ulation  of  Eqs. (41)and(42)  by  rederiving  it  for  the  case  where  the  two  primary  waves 


are  at  an  angle  relative  to  each  other.  In  this  Instance  a  resonant  interaction 
should  be  possible  at  least  in  parts  of  the  field  at  the  frequencies  of  interest 
(note  that  the  interaction  is  considerably  more  complicated  than  that  due  to 
simple  plane  wave  fields)  but  less  and  less  resonant  at  other  frequencies  in  the 
nonlinearly  generated  spectrum. 


iaaiiiaiiaiaai* 
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3.  Results  and  Statement  of  Continuing  Research 

Before  proceeding  to  a  computer  investigation  of  the  mathematical  models 

outlined  in  this  report  it  was  felt  desirable  to  obtain  upper  bound  estimates 

of  the  extent  to  which  the  conversion  efficiency  of  a  parametric  array  can  be 

Increased  when  all  frequency  components  other  than  the  primary  waves  and  the 

difference-frequency  are  suppressed.  Making  use  of  a  computer  program  de- 

49 

veloped  by  McKendree  and  Fenlon  which  operates  in  the  time  domain  but  is 

subject  to  Iterative  spectrum  analysis  this  comparison  is  shown  in  Figs.  2 

and  3  for  the  case  of  a  virtually  lossless  plane  wave  parametric  interaction 

in  which  the  frequency  ratio  of  the  primary  waves  is  3/2.  In  this  instance, 

it  can  be  seen  from  inspection  that  the  three  frequency  case  has  a  maximum 

difference-frequency  amplitude  nearly  an  order-of-magnitude  than  that  of  the 

all-frequency  component  spectrum.  Moreover,  the  first  peak  of  the  difference- 

frequency  in  Fig.  2  occurs  at  o'  »  3.5  at  R  »  3.5/aQ.  Hence  at  a  value  of 

0Q  between  j,  and  10  it  can  be  seen  from  Fig.  4,  which  was  first  evaluated  by 

47 

Fenlon  and  McKendree,  that  the  beamwidth  is  almost  fully  established.  This 
ensures  that  if  the  medium  in  which  the  three  frequency  interaction  is  term¬ 
inated  an  ideal  and  much  more  efficient  parametric  array  will  be  formed  rel¬ 
ative  to  that  which  would  have  existed  under  normal  conditions.  Similar  but 
even  more  efficient  comparisons  are  then  shown  in  Figs.  5  and  6,  for  »  6/5 

and  in  Figs.  7  and  8  o^/u^  -  11/10  respectively.  These  cases  provide  typical 
upper  bounds  to  the  maximum  realizable  gains  in  conversion  efficiency  achievable 
via  absorption  and  dispersion  mechanisms,  both  physical  and  artificial. 

The  the  present  time  we  are  in  the  process  of  programming  and  testing 
frequency-domain  computer  methods  to  determine  the  effect  of  dispersive  mecha¬ 
nisms  in  layered  media,  and  using  then  to  obtain  the  on-axis  solutions  of 
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f inice  beams  in  relaxing  fluids.  The  latcer  is  particularly  difficult;  to 

solve  numerically  because  the  diffraction  losses  are  extremely  difficult  to 
include.  These  are  however,  defined  paraxially  by  Eqs.  (41)  and  (42) 

were  derived  within  and  therefore,  further  analysis  and  simplification  of 
these  results  is  expected  to  reveal  the  correct  way  to  include  them  in  a  step 
by  step  wave  propagation  analysis. 

Once  the  computer  models  are  running, sets  of  nnndlmensional  plots  of  the 
parametric  fields  ,in  dispersive  media  can  then  be  evaluated,  plotted  and  stud¬ 
ied  in  detail  for  subsequent  relation  to  real  and  artificial  acoustic  media. 

Before  concluding  this  discussion  it  is  interesting  to  check  the  accuracy 
of  Figs.  2,  5,  and  7  by  comparison  with  the  exact  three  frequency  solutions 
given  in  Appendix  A.  As  far  as  the  first  peaks  of  the  latter  are  concerned 
they  are  in  excellent  agreement  with  the  analogous  at  exact  curves  of  Figs. 

(Al) ,  (A2),  and  (A3).  Subsequent  deterioration  of  the  computed  curves  Id  due 
both  to  the  presence  of  small  absorption  losses  (i.e.  -  10^)  and  to  "nu¬ 

merical  error".  However,  in  general  only  the  first  peak  of  the  difference- 
frequency  is  of  any  interest  because  at  this  point  the  dispersive  medium  must 
be  terminated  and  the  amplified  difference-frequency  signal  released  into  the 
unbounded  surrounding  medium.  Hence  the  utility  and  validity  of  the  numer¬ 
ical  methods. 

During  the  next  phase  of  the  contract  as  well  as  performing  the  numerical 

analysis  referred  to  above  we  will  carry  out  an  analytical  investigation  of  finite- 

amplitude  waves  of  Gaussian  cross-section  intersecting  at  a  f  '.xed  angle  in  dispersive 
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fluids.  Unlike  the  simple  plane  wave  case  Investigated  by  Rudenko  and  Solugan 
we  do  not  know  if  this  will  result  in  a  resonant  interaction  at  all  points  of 
the  field  or  only  in  the  far-fleld  of  the  interaction.  This  is  one  objective  of  the 
investigation.  The  other  is  to  establish  how  rapidly  remnant  interaction  falls 
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of£  throughout  the  field  at  nonlinear ly  generated  frequencies  other  than  the 
difference-frequency,  thus  determining  the  extent  to  which  the  latter  can  grow 
due  to  dispersive  coupling. 


3/2  BIFREQUENCY 


figure 
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Appendix  A 

If  and  ?2  are  plane  primary  wave  amplitudes  of  frequencies  and 
respectively  normalized  with  respect  to  P10  and  P„  is  the  normalized  differ¬ 
ence-frequency  (i.e.  ui_  *  u>^  -  uij)  field  resulting  from  their  nonlinear  inter¬ 
action  then  in  the  absence  of  lossos  and  other  frequency  components  if  the 
three  waves  have  a  phase  of  7r/2  then 


— i.  p  p 

2  *2  e- 


1N, 

— —  p  p 

2  *r~ 


where 


dP_  IN.  „ 

d?-  ’  -  T  P1PZ 


(*>,  U), 

and  N1  -  ^  ,  N2  -  uj  ,  N.  -  ^  -  Nj.  (A5) 

The  normalized  range  R  *  ,  where  z  is  the  half-Rayleigh  distance  of  the  projec- 

1  0 
o 

„  2 

z  p0  co 

tor  and  a  ■  —  ,  where  z  •  w= — r  is  the  critical  range  for  a  mean  wave- 

0  C  P  r « 

^  c  10  o 

number  kQ  *  y  (k^+kj) .  As  shown  by  Bloembergen,  EqB.  (Al)  to  (A3)  have  an 
exact  analytical  solution  in  terms  of  Jacobian  Elliptic  functions,  summarized 
by  Scott^. 

|PJ  -  ?1kj  cdCkjO' ,  k2)  (A6) 

|P2I  *  1*20  nd(kia'*  (A7) 


P"  ”  lN2  I  ^  P2°  k25 
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